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Abstract. The Grassmann structure of the critical XXZ spin chain is studied in 
the hmit to conformal field theory. A new description of Virasoro Vernia modules 
is proposed in terms of Zamolodchikov's integrals of motion and two families of 
fermionic creation operators. The exact relation to the usual Virasoro description 
is found up to level 6. 



1. Introduction 

In the present paper we continue the series of works [H |2l |3] on the XXZ model. In 
[5] we considered it in the presence of the Matsubara direction, or equivalently the 
six vertex model on a cylinder. We computed the normalised partition function with 
a defect localised between two horizontal lines, which corresponds to an insertion of 
a quasi-local operator: 



TrsTrM(Ts,Mg2-^+2«5(o)j 



Here q = e^^^ is related to the coupling parameter (see (12. 2p below), and T^^m stands 
for the monodromy matrix on the two tensor products of evaluation representations 
of Uq{s[2)'- one for the horizontal (or 'space') direction S, and another for the vertical 
(or Matsubara) direction M. For more details see section [2] below, in particular fig. 
1. It was important in [3] to incorporate inhomogeneities in the Matsubara chain. 
This allows, for example, to consider the temperature expectation values in the spirit 
of [U [5], by adjusting inhomogeneities and taking the limit to the infinite chain in 
the Matsubara direction. 

The clue to our calculation was the introduction of operators t*(^), b*((^), c*(^) [2J 
which, by acting on the primary field g^^'^C^) , create the space of quasi- local operators 
on the horizontal chain. More precisely, quasi-local operators are created by Taylor 
coefficients of t*(C), b*((^), c*{Q at the point C,"^ = 1. In this paper we change 
the definition of b*((^), c*(^) from those of [21 [3] by applying certain Bogolubov 
transformation. Compared with the original ones, they have better asymptotic 
properties. We shall explain this in section |2J 

There is an obvious similarity with conformal field theory (CFT), where the de- 
scendants are created from the primary field by the action of the Virasoro algebra. 
Our aim in this paper is to examine the scaling limit of our construction in the 
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critical regime, and to establish its precise relation with CFT. We shall consider the 
case of a homogeneous Matsubara chain. 

The functional fll.ip is non-trivial only for operators of spin zero. As it turns 
out, for the study of the scaling limit, it is quite useful to relax this restriction. We 
shall first introduce the following generalisation of fll.ip which is of interest on its 
own right. For s > we define 

, , ,,,, , TrsTrM(l4-^^Ts,M q'^' h*^^,_, ■ ■ ■hl^^.iq'-'^^O) 



TrsTrM(ri-^^Ts,M q^^' b:,,,_i ■ • ■ b;,,o(g2"5(o)) 

where the .^'s denote the coefficients of the singular part of b*(C) at = 0. 
When s < 0, a similar definition is in force using the expansion coefficients of c*{Q. 
As long as is taken generically, this definition is independent of its choice (see 

section [2] for more details). In general the dependence on enters, but only in 

a "topological" way. 

Needless to say, what we we are dealing with is a lattice analogue of the screening 
operators a la Feigin-Fuchs-Dotsenko-Fateev [6l [7] . It is interesting to see that quasi- 
local operators and screening operators both arise from the same operators b*(C), 
c*(C), as expansions either around (^^ = 1 or (^^ = 0. 

After these modifications the main formula of [3] remains valid. It reads 

(1.2) z-^^{t*{o ■ ■ ■t*(c;)b*(c+) ■ ■ ■b*(C)c*(C7) ■ • ■c*(cr)(9'"^^°^)} 

p 

where the functions p{Ci\i^^ n + a,s) and oJ^Ct ^ ^) ^'^^ defined by the data of 

the Matsubara direction. We refer to this formula as the determinant formula. Due 
to the Bogolubov transformation of b* and c* the function a;(C, a, 0) is slightly 
different from a) used in |31 IB]- 

Now let us turn to the scaling limit. We have two twisted transfer matrices 
TmiCi 1^ + ot) and Tm(C5 '^) in the Matsubara direction. As the the number of sites 
n becomes large, their Bethe roots tend to distribute densely on M+. Fixing i? > 
and introducing the step of the lattice a, we consider the limit 

(1.3) n — 7- oo, a — 7- 0, na = 2'kR fixed . 
At the same time we rescale the spectral parameter as 

(1.4) C = {Cay\ A fixed, 

so that the Bethe roots close to stay finite in terms of the variable A. Here C 
is a constant chosen for fine tuning (see section |8l (18. 3p ). In this limit the twisted 
transfer matrices turn into the transfer matrices of chiral CFT on the cylinder Cyl = 
C/27rii?Z, introduced and studied by Bazhanov, Lukyanov and Zamolodchikov (Qj 
[To] . We wish to mention here that the present work owes a great deal to these 
remarkable papers without which it would have been impossible. Details about the 
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scaling limit can be found in section |H] below. The relevant CFT has central charge 

c = 1 



1 - V 

We shall parametrise the conformal dimension as 



Z/2 



.((«- 1)2-1), 



" 4(1 -z.) 

and write the action of the Virasoro algebra on a local field <^{y) as 

In the following we set a = Ca, and let limgcaiing indicate the scaling limit fll.3p . 
(I1.4p . The functions entering the determinant formula (II ■2p also have finite limits, 

P^'^{X\k, k') = lim p(Aa'^|K, a, s) , 

scaling 

4 k' ,a) = lim u{Xa'^ , fia'^\K,a, s) , 

scaling 

where 

(1.5) K=K + a + 2^s. 

So all these partition functions (11. 2p have finite limits. They should have some 
definite meaning in the context of CFT. We contend that they are the three point 
functions of the descendants of the chiral primary field 0a(O), computed in the 
presence of two other primary fields (or their descendants) inserted at the two ends 
of the cylinder. 

More specifically, we conjecture that the following picture holds true. First, the 
creation operators tend to a limit, 

2t*(A) = lim t*(Aa^), 2/3* (A) = lim h*{Xa''), 27*(A) = lim c*{Xa'') . 

scaling scaling scaling 

As A — OO, these operators have asymptotic expansions of the form 

oo 

(1.6) log(r*(A))^^T;^._,A-^, 



oo 



T*(A) ^ Vt*(A) ^ 



In the limit, the quasi- local operator g2Q!5(o) ^gpomes the product of two chiral 
primary fields 0o(O) (8>0„q,(O). The operators T^j^i and the quadratic combinations 
/32j_i72fc-i S'Ct only on the left component 0q(O), and create the entire Verma module 
spanned by the Virasoro descendants 

l_mi ■ ■ ■l-m,(0a)(O) . 

Furthermore, if Y-^^^ is chosen to be generic then 

(1.7) hm z-'^{t*iO ■ ■ •t*(c;)b*(c+) • • ■b*(C)c*(c) ■ ■ •c*(cr)(g'"^(°))} 

scaling 

= 2^+2^zr'{T*(A?) ■ • • r*(A;)/3*(A+) ■ ■ • /3*(A+)7*(A;) • ■ ■ 7*(Ar) (0.(0)) } • 
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In the right hand side, the symbol Z'^'^ {X(0)} stands for the three point function 
normahsed as Z'^^ {0a(O)} = 1, with X(0) inserted at a; = and the primary fields 
being inserted at x = oo and x = — oo, respectively (see fl6.8p below). 

i—s) 

Non-generic choice of Y-^ corresponds to replacing the primary fields at x = ±oo 
by their descendants. In this paper we discuss only the case of generic Y^^\ 
The coefficients in (11.61) are homogeneous operators in the sense that 

[ lo , t;,_i] = (2j - l)r;^._i, [ lo , /3;._i7;,-i] = {2i + 2j - 2)l3l,_,^l^_, . 

Hence, for each degree, the descendants created by I_fc's, and those created by 
T2j_i's, /32j-i's and 72j_i's, must be finite linear combinations of each other. The 
main goal of this paper is to show, for low degrees, that this is indeed the case, and 
that the coefficients can be found explicitly. 

To determine the coefficients of the linear combination, we compare the values 
of Z^'^ . For the Virasoro descendants, they can be easily computed by the con- 
formal Ward-Takahashi identities. For the descendants by T2j_i and others, we 
need the coefficients of the asymptotic expansion of the functions p/j(A|/t, k') and 
ujr{\, k', a) . In section[TD]we develop a systematic method for computing them. 
We note that in both cases the results are polynomials in the conformal dimensions 
A^+i, Ak/_|_i. We may regard them as independent variables and compare the coef- 
ficients, since s in (II. 5p can take any integer values. This was one of reasons for us 
to introduce the screening operators. 

We consider first Tj^.i. CFT allows an integrable structure based on Zamolod- 
chikov's integrals of motion i2m,-i |Il]- With the above procedure we are led to a 
result which should not be surprising, 

''"2m- 1 ~ ' 12771- 1 ) 

where Cm are some i^-dependent constants which can be found in [lOj . 

We then consider the action of /32j_i's and 72j_i's. Because of a technical difficulty 
we have not been able to compute the asymptotics of ur{X, fi\K,, K,',a) for k ^ k! . 
Here we restrict to the case n = n'. Since Z^''^(i2„_i(X)) = for any X, restricting 
to K = k' means that we consider the quotient space of the Verma module modulo 
the action of the integrals of motion. We assume that the vectors 

i2fci-l ■ ■ ■ i2fcr-ll-2mi ■ ■ ■ l-2ms (0Q:(O)) 

span the Verma module, so the quotient space is created by the l-2m's. With primary 
fields as asymptotical states, we can compare up to the level 6. It should be added, 
however, that up to this level the system of equations is overdetermined. So the very 
possibility of finding a solution is the strongest support of our fermionic picture. 
We give one example on the level 4: 

PllliMO)) = \D,{a)D,{2 - a) {\\ + ^-^i^u) (0„(O)) , 
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where 

da = |V(25-c)(24A„ + l-c), 

In the right hand side, we have a particular combination of the Virasoro descendants. 
This equation says that its three-point function remains of the same determinant 
form before and after integrable perturbation. 

The text is organised as follows. In section |5] we review the results of [3] and 
describe the Bogolubov transformation mentioned above. In section [3] we define the 
functions p and u. In section H] we define screening operators on the lattice, and 
describe a generalisation of the previous results. In section |5] we start discussing 
the scaling limit of the XXZ chain, examining the behaviour of the Bethe roots in 
the Matsubara direction as the length of the chain becomes infinite. Sections [6] and 
[7] are a review of the CFT integrals of motion on the cylinder, and the series of 
works of Bazhanov, Lukyanov and Zamolodchikov (BLZ). We explain in section |8] 
how the Matsubara transfer matrix turns into that of BLZ in the continuous limit. 
In section IH] we discuss the CFT interpretation of the scaling limit in the space 
direction. In section [10] we study the asymptotics of Thermodynamic Bethe Ansatz 
(TBA) function a for CFT. In section [11] we find the asymptotical expansion of uj 
for K = k'. In section [12] we compare descendants created by l-2m with those created 
by (3*2j_iS and 72j_i and give some concluding remarks. In appendix we present 
general properties of asymptotics of u which apply to the case k ^ k,'. 

2. Review of previous results 

Let us start with a brief review of the papers [H [2], [3] . Consider the XXZ spin 
chain in the infinite volume. The space of states of the model is 

oo 

Sjs= (g) , 

j=-oo 

and the Hamiltonian is given by 

oo 

(2.1) H=l J: {alal^, + alal^, + Aalal^,) , A = + q~') . 

fc=— oo 

We consider the critical XXZ model in the following range of the coupling, 

(2.2) q = e"^^ l/2<u<l. 

Together with S^s we consider the Matsubara space i^M- In [2] the most gen- 
eral case was treated: namely, i^M was the tensor product of spaces of different 
dimensions, and to every site m an independent inhomogeneity parameter was 
attached. In the present paper we shall restrict ourselves to the case 

n 
j=l 
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We consider the monodromy matrix Ts,m- Mathematically this is nothing but 
the image of the universal R matrix of Uq{5l2) on the tensor product of evaluation 
representations corresponding to Sjs and S)m- It has been said that we shall consider 
a homogeneous Matsubara chain only. In the notation of [3], this correspond to 
setting Tm = q'2 for all m. We shall absorb this into redefinition of the L-operator 
comparing to |3j. Let us write the definition explicitly: 

r\ 

oo 



where 

n 
m=l 

with 

LUC) = Q"'"^"- - CV^"^"- - Ciq - + o-ot:) . 

A local operator on i^g is by definition an operator which acts nontrivially only 
on a finite number of the tensor components of i^s- More generally we consider 
quasi- local operators with tail a, which are operators of the form 

j=-oo 

with being local. In this notation S = S{oo) is the total spin. In [3] we computed 
the expectation values defined by 

TrsTrM(Ts,Mg"^^+'°^^°)0) 
(2.3) Z'^ 1^2-^^(0)0 1 = \ ^ . 

^ ^ TrsTrM(Ts,Mg2'^^+2"^(°)) 

This is a linear functional on the space Wq^q of spinless quasi-local operators with 
tail a. 

It is helpful to think of the functional Z'^ as a ratio of partition functions of the six 
vertex model on the cylinder. For example, the numerator of (12.31) can be presented 
graphically as follows. 
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u 

b 
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r 
a 



= Lij 



x=q' 



Fig.l: Partition function on the cylinder. 

The functional Z'^ is a ratio of two partition functions on the cylinder. On each 
crossing of a row and a column, one associates the Boltzmann weights of the six 

3 

vertex model. On a particular row there are also twist fields q'^^^ ( marked by crosses ) 
and q^°'^'^'>"j (marked by circles). The numerator of corresponds to a lattice with 
defects representing an insertion of a local operator. 



On this picture the summation is performed over all edges except the broken ones 
in the middle. The arrows on the broken edges are fixed, representing the local 
operator 0. The one dimensional sublattice going in the infinite space direction 
will be called the space chain, while the compact one dimensional sublattice in the 
Matsubara direction will be referred to as the Matsubara chain. 

We computed the expectation values (12. 3p using the fermionic description of Wq^q 
found in [2]. Let us briefiy recall it. Consider the space 



S,S 1 



where 'Wa~s,s denotes the space of quasi-local operators of spin s with tail a — s. 

We have the creation operators t*((^), b*(C), c*{() and the annihilation operators 
b(C), c{() which act on W*^"). To be precise the operators t*(^), b*(^), c*(() were 
defined in [2] as formal power series in (^^ — 1, the quasi- local operators in W*^") are 
created by coefficients of these series. However, when the series are substituted into 
Z'^ the result allows analytical continuation. So, in the present paper we shall adopt 
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another point of view which is similar to that of CFT. Namely, we shall consider 
the operators t*((^), b*((^), c*{() as analytical functions. Then the relation to real 
quasi- local operators is achieved by considering t*(^), b*(C), c*(^) around the point 

The creation-annihilation operators have the block structures 



' a—s,. 



(2.4) t*(C) : W„„,,, ^ W, 

b*(C),c(C) : W^-s+l,s-l ^ V^a-s,s, 

c*(C),b(C) : W„_,_i,,+i ^ W,_,,,. 

The operator r = t^/2 plays a special role: it is the right shift by one site along the 
space chain. 
We have 

b(C)(g2"^(°))=0, c(C)(g'"^(°)) = 0, 

[c(o, c*(c)] + = ^(e/c, «), [b(o, b*(c)] ^ = -^(c/e, «) , 

where 

(2.5) ^{C,a) = C: + ^ 



2{e - 1) 

The operators in the space Wq, o are created from the primary field by action 

of t*'s and of equal number of b*'s and c*'s. The completeness [I2] says that the 
entire space Wo,o is generated by coefficients of the creation operators considered as 
series in — 1. Certainly, this description is reminiscent of CFT. 

The main result of [3] is the following relations which allow for recursive compu- 
tations of the expectation values: 

(2.6) Z^{t*iC)iX)} = 2piC\n, K + a)Z^{X} , 

(2.7) Z^{h\aX)] = uj{Ci\n,a)Z^{c{aX)]^ . 

r 

(2.8) Z«{c*(C)(X)} = -^ya;(e,CI«:,«)Z'^{b(0(X)}^, 

r 

where F goes around all the singularities of the integrand except = C^- We think 
no further explanation is needed nowadays when the method of CFT is a part of 
common knowledge. 

The functions p(C|k, k + a) and c<;(C,^|/t, a) will be defined in section [3l We 
changed the notation for the former from to [3]. The present notation agrees better 
with the explicit formula given below. The set of equations (12. 6p . (12. 7p . (12. 8p implies 
a determinant representation for the expectation values ^, we shall discuss this 
later. 

Let us describe the modification of b*((^), c*{Q by a Bogolubov transformation 
which was mentioned in the introduction. Denoting the operators used in |3] by 
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^rat (C), c:^t(C), define 

(2.9) b*(C) = b:,,(C) + / Dc^5A-V(C/e,«) ■ c{0% , 

r 

c*(C)=<,(C)-^ I /^c^^A-V(e/C,«)-b(0^. 

r 

where D^^ is the following finite difference operator of the second order, 

^c/(C) = /(Cg) + /(Cg-')-t*(C)/(C). 

The function A'^'^ip{(, a) is transcendental. For 

> 0, --<Rea<0 



we define it as 

oo 

(2-10) ^-V(C,«) = -VP I ' a^ ^fa")o^.-.. 



'(l + (C/r/)^) ''""'27rir/2 



where the principal value is taken with regards to the pole at = 1. In general 
we define it by analytic continuation with respect to both a and C^, obtaining a 
meromorphic function of logC- It is bounded at logC — ?■ ±oo, and its singularities 
closest to the real axis are the simple poles at log ( = ivriz/ with residues of opposite 
signs. 

The function D(^D^A'^'^^{(/$^, a) is regular at C = ^, so the Taylor series for b*((^), 
c*{() at = 1 are well-defined. The function u changes following the change of b*, 
c*: 

(2.11) u{C, a) = ^rat(C, el^t, a) + D^D^A^^^iC/^, a) , 
where 

(2.12) DJiC) = fiCq) + fiCq-') - MC\f^, + «)/(C) • 

Let us mention one marvellous property of our modified operators b*, c*. 
The main subject of our original study [H |2] was the following normalised matrix 
element 

(2.13) z.{,-.».0}^ip2«^, 

(vac I g"'"''^^-' I vac) 

where |vac) denotes the ground state of the XXZ Hamiltonian in the infinite volume. 
Equivalently, the numerator of fl2.13p is the partition function of the six vertex model 
on the plane with a defect localised between two horizontal lines. 

It is easy to see using the formulae from [H [2] that, if we create the quasi- local 
fields by the operators t*{(), h*{(), c*(^), then Z^o vanishes on all of them with the 
sole exception of the descendants created by r = tl/2. For the latter we have 

^oo{T'^(g^°'^^°^0)} = 1, m G Z . 
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We know from the algebraic construction [2j that the Taylor coefficients of the 
part b*g^^((^), c*g^^(^) in (12.91) produce only rational functions of q and q". Hence all 
transcendental pieces of the expectation values (12. 9p come from rewriting a given 
quasi-local operator in the fermionic basis and picking Taylor coefficients of the 
function D^D5A^V(C/^, «)• 

The property of t*{Q (except t^) and b*(^), c*{Q mentioned above is analogous to 
that of the Virasoro generators in CFT: on the Riemann sphere, all normalised one 
point functions of the descendants vanish due to the conformal invariance. Strictly 
speaking, the one point function of the primary field also vanishes, so a word of 
clarification is necessary. Take a massive model with a mass scale m and consider 
the conformal limit m — )■ 0. While the normalised one point function of the primary 
field stays equal to 1, those of the descendants vanish in the limit, because for 
dimensional reasons they contain additional powers of m. 

At this point one may wonder why we did not do the Bogolubov transformation 
killing completely the function u in (12.71) . (12.81) . The answer is that it is impossible 
to rewrite the left hand sides of (12. 7p . (12. 8p because the function u cannot be written 
as a function of p. 



3. Functions p and uj 
Introduce the twisted Matsubara transfer matrix: 

(3.1) TM(C,«:) = Tr,(r,,M(C)g'"^^'). 

Let |k) be the eigenvector of T{1^k) whose eigenvalue is maximal in the absolute 
value. Similarly let (ft + a\ be the eigencovector of T(l, k + a) whose eigenvalue is 
maximal in the absolute value. It is well-known that these eigenvectors have spin 0. 
We call them the maximal eigenvectors, and assume that they are not orthogonal. 
We denote the eigenvalues of Tm(C? i^) (resp. Tm(C) ^^ + Q^)) oii (resp. {k + a\) 
byT(C,K) (resp. T(C,K + a)). 

Then p(C|'^, k + a) is defined by T((^, k) and T(C, k -|- a). We have 

(3.2) p(C|fi;,K + a)- 



T(C,«:) 

The function uj{C,,^\K,a) is more complicated. In [5] it was shown that it is com- 
pletely determined by two requirement: the singular part and the normalisation 
condition. Then they were explicitly solved in terms g-deformed Abelian integrals. 
In the present paper it is convenient to use an alternative, TBA-like, description 
used in [8]. Let us explain this simplifying a little the notation of [8]. 

Together with the transfer matrix Tm(C) ^) we consider in [3] Baxter's Q-operators 
Q^(C,k). In this paper we use only one of them: Q^{C,ii^) denoting it just as 
Qm.{.CtI^)- It is defined by the trace over the highest weight representation of the 
g-oscillator algebra with generators a, a*, D (see P] for notation). 



(3.3) Qm{C^) = C''^'^il-q-''^''-'^^)Tr-{To,,^^{Oq-'-''-) 
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n 



^Osc,m(C) — Y\. -^Osc,m(C) 



m=l 

Its eigenvalue on the eigenvector discussed above will be denoted by Q{(,k). The 
main role in TBA is played by the function: 

a(C)Q(Cg 

where 

«(c) = (i-gcr, dic) = ii-q-'er- 

It follows from the Baxter equation that the solutions to the equation 

a(C,«^) = -1, 

are the zeros of either Q{C, k) or T{(, k). The function a(C, f^) satisfies the nonlinear 
integral equation [131 IH E] ■ 

(3.5) log o(C, k) = -2mvK + log - / log (1 + «)) 



7 



where the cycle 7 goes around the zeros of Q{C, n) (Bethe roots) in the clockwise 
direction, as opposed to all other contours. 

We shall need slightly more general kernel than K{(/C,), so, let us define them 
together. First, we introduce operations: 

(3.6) Ac/(C) = fiCq) - fiCq-') , 

5c"/(0 = /(Cg)-p(CI'^,'^ + «)/(C). 

Then 

(3.7) ir(C,a) = -^Ac^(C,a), K{C) = K{C,0) . 

Ztxi 

We shall use the the following notation: 

1-^9= f{v)9{v)dm{r]) , 



where the measure is given by 
(3.8) dm{r]) 



rfp{rj\K, K + a) (1 + a(?7, k)) ' 
Now we introduce the resolvent of certain integral operator 

(3.9) -Rdress -^dress ^ ; 
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where stands for the integral operator with the kernel K{(/C,,<y). Introducing 
further the two kernels 

(3.10) /ieft(C,0 = Li '^c"V^(C/e,«), /nght(C,0 = 5^i^{C/^,cy) , 

we are ready to write the definition of [8] cleaning it from irrelevant auxiliary objects 
and taking into account the modification (12. lip : 

(3.11) icu(C, a) = (/left -k /right + /left * ^drcss * /right) (C, 

+ 5-5-A-V(C/e,a). 

4. Introducing screening operators on the lattice 

Now we want to describe an important generalisation of results of [3J. We have 
three constants: a which defines the tail of the operator and k, k + a which define 
the twist of Matsubara transfer matrices at +00 and —00. As it has been said in the 
introduction we need more freedom. Let us explain how an additional parameter 
s G Z can be introduced into Z'^. We shall see later that in the scaling limit 
introduction of this parameter leads to emancipation of k + a from k and a. 

Let us consider the trace 

(4.1) TrsTrM(vi-^)rs,Mg^''"+^("-^)"(°)0(^)) , 

(—s) 

where Y-^ carries spin —s. For definiteness we suppose s > 0. It follows from the 
ice condition that in this situation the operator 0*^*^ must have spin s. 

We assume that among the eigenvectors of the transfer matrices Tm{C, i^) and 
^m(C) i^+o — s) there are maximal ones \k), \K+a — s, s) with eigenvalues T{(, n) and 
T[(, K + a — s,s) which are defined by the requirement that T(l, k) ■T(1, n + a~s,s) 
is of maximal absolute value among all the pairs of eigenvectors. We assume further 
the generality condition: 

(4.2) {K\Y^'^\K + a- s,s) ^0 

Obviously the difference between spins of \n + a — s, s) and |k) must be equal to s. 
We make the technical assumption that spin of \k) remains equal to zero. 

The natural idea is to create the operators g2(a-s)s(o)Q(s) c*, t* having an 

excess of operators b*: 

g2("-^)^(o)o(^) = b*(G)---b*(^,) 

X b*(C) ■ • ■b*(C)c*(cr) • ■ ■c*(C)r(c?) ■ ■■t*{0{q'''"^'^) . 

However, the formulae (12.71) are not applicable in this case. Let us explain why it 
is so. 

The method of [3] requires to start the consideration by the operator b*(^i) which 
is the closest to Ts^mi'^'^^ ■ The operator ^^'^h*^^{^i){X) is a meromorphic function of 

with singularities at the points (Cj~)^- It satisfies n+1 normalisation conditions [3]. 
The additional term in (12. 9p is of the form satisfying (12. 7p from the very beginning. 
These singularities and normalisation conditions are studied algebraically, they do 
not change comparing to [3] where the spin of X was equal to —1. However, the 



GRASSMANN STRUCTURE IN XXZ MODEL 



13 



behaviour at zero changes: in the present case the spin of X equals s — 1, and 
according to [2] ^'^"'h*{^i){X) does not vanish at zero. Generally, 

s-l 

(4.3) C""b*(C)(X) = Y,C''h*^M) + C""b;,g(C)(X), X e W,_,+i,,_i , 

j=0 

where C~"b*gg(^)(X) vanishes at zero. But in [3J the fact that C""b*(C)(X) vanishes 
at zero for spin(X) = —1 was important when deriving (12.71) . As a result in the 
present case we do not have enough conditions to define u. 

Let us turn this problem into advantage. The operators b^^, j = 0, ■ ■ - s — 1 
constitute a finite Grassmann algebra. So, we just consider — )■ and replace 
b*(^i) ■ • • b*(,^s) by bj^ s_i ■ ■ ■ b^ q. Now move one of the remaining operators b*, 
namely, b*(C;^'") to the left. Obviously, (C]^)~"b*(C]^) vanishes as (Ci")^ because 
the singular part disappears due to multiplication by bj^ s-i • • ■ bj^ q. Effectively, 
this operator reduces to b*gg(^]^). 

It is not hard to see that we can introduce a;(C,C) and obtain (12. 6p . (12. 7p . (12. 8 p 
for the functional on Wq, q: 

, . TrsTrMflir^^Ts.M g^'^^b^oo.s-i ■ ■ -h* ^^^{q^^^^^^O) 

The function p changes in the most natural way to 

T(^, K + a — s, s) 



(4.4) p(C|/€,/€ + «,s) 

The function c<;(^, a, s) is defined by ( 13. lip , replacing p(C|'^, f^ + o) by p(C|/t, + 
a, s) but keeping the same a(C, «:). 

Let us discuss one important property of Bethe vector of spin s. The basic object 
in the theory are the transfer matrix Tm(C5^ + a — s) and Baxter's Q-operator 
Qm(C; k + a — s). Their eigenvalues on the vector \k + a — s^s) have the following 
analytical properties: T(C, k + a — s, s) is a polynomial of of degree n, 

g(^, ^ + a-s,s)= C-"-'^+'M(C, K + a-s,s), 

where A{C,, n + a — s,s) is a polynomial in (^^ of degree n/2 — s. In terms of T and 
A the Baxter equation reads 

(4.5) T(^, K + a — s, s)A{C, k + a — s, s) 

= g-"'ci(C)v4(Cg, K + a-s,s) + q^'a{C)A{Cq~\ K + a-s,s). 

with 

(4.6) k' = K + a + 2^s. 

Formally, k' is defined modulo TLjv. However, the eigenvalues are multi- valued 
functions of k' and we have to be careful about the choice of the branch. The choice 
in (14. 6 p is consistent with the semi-classical limit i/ — )■ 1. We shall return to this 
point when we discuss the scaling limit. 
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Notice that if we write the Baxter equation for twist and spin 0: 

(4.7) r(C, k')A{C, k') = q-'-'diOAiCq, k') + q''' a{OA{Cq-\ k') , 

it looks exactly the same as (14.51) if we identify T{(, k') with T{(, k + a — s, s) and 
A[<!^,k') with A{(,K + a — s,s). There is, however an important difference: the 
polynomial A{(, k') in ( 14. 7p is of degree n/2 while the polynomial A{(, k, + a — s,s) 
in (14.51) is of degree n/2 — s. Still, this similarity will be very important for us later 
when we shall discuss the scaling limit. 

Similarly to the previous discussion we can consider an operator Y-^^ which carries 
positive spin s. Then the operator c*{Q will have nontrivial behaviour at ("^ = 0: 

s-l 

(4.8) Cc*(C)(x) = Y,c''c*^AX) + c"<g(c)(x), X e w„+,_i,_,+i , 

j=0 

and we can repeat the entire procedure using c^^j instead of b^^. So, s in (14. 6 p 
can take any integer value. 

Obviously, what we are doing here is nothing else but introducing screening op- 
erators on the lattice. This is important for relating to the CFT. The screening 
operators anticommute among themselves, the same is true for C^j- Naively, 
one could say that b^j anticommute with j, but this does not make sense be- 
cause the product of these operators do not act nontrivially on any subspace of 

Remark 4.1. Let us mention one more, less dramatic, generalisation of the results 
of |3j . Clearly the functional Z'^'''^ is independent of the choice of provided the 
condition (14. 2 p is satisfied. However different choices are also possible. For instance 
one can take any eigenvector \A) and eigencovector {B\ of the Matsubara transfer 
matrices and consider the projector \A){B\. The main formula is applicable in this 
more general setting. 

5. Scaling limit 

We take a, k to be real, and restrict our consideration to the region 

(5.1) \k\ < 1, \k'\ < 1. 

Then we continue analytically. 

The crucial data for our construction are the eigenvectors \k + a — s,s), \k,). 
Consider the second of them. The corresponding eigenvalue T{(, k) corresponding 
to maximal eigenvector |k) of spin zero. Denote by Q{(,k) the eigenvalue of the 
Q-operator on \k). These eigenvalues have the form 

n/2 

Q{c,^) = c''ii(^-c/e,), 

p=i 

n 

nc,^) = {q^+q~nU(^-eK)- 
p=i 



GRASSMANN STRUCTURE IN XXZ MODEL 



15 



In the domain \k,\ < 1, the vector |k) is uniquely characterised by the two require- 
ments for the roots: G M, and > > 6'| for all k,j. Let us study the 
behaviour of the Bethe roots as n — )■ oo. Suppose they are numbered in the order 
^i<^2<^3<'''- the limit n — ?■ oo they are subject to the Lieb distribution 
[T5l : for 1 ^ m we have 



(5.2) ^-f^-^ = "iid + fJ)+0{j-, 

We are interested in the Bethe roots which are not very far from ("^ = 0. In other 

words we assume 1 ^ m ^ n. Since is small, we can drop the term in the 
(15.21) . obtaining 

(5.3) ^ (vr-j . 

Similar power law is obeyed by 9m- 

So far we have concentrated on the ground states in Matsubara direction. But 
according to Remark 14.11 the main formulae can be generalised to arbitrary Bethe 
states. There are low- lying excited states which satisfy (15. 3p . and to which the same 
analysis as for the ground states apply. Readers who are famihar with the papers 
by Bazhanov, Lukyanov and Zamolodchikov [9], [101 [IB] would immediately say that 
in the limit we obtain the CFT transfer matrices treated by them. We shall come 
to this relation later in section [71 

The formulae (12. 6p . (12. 7p and (12.80 imply the explicit expression 

(5.4) z-'^{t*iO ■ ■ ■t*(c;)b*(c+) ■ ■ ■b*(C)c*(c-) ■ ■ ■c*(cr)(g'"^(°^)} 

p 

= Y['^p{Ci\i^,f^ + a,s) xdet (w(C+,C7l'«,tt,s)). .^^^^ . 

i=l 

We consider the scaling limit in the Matsubara direction, 

(5.5) n — i- oo, a — )■ 0, na = 2ttR fixed . 
Let us introduce the following strangely looking notation 

(5.6) a = Ca , 

where C is some i/-dependent constant which will be needed for fine tuning compar- 
ing the scaling limit with CFT. 

The following limits exist for finite A: 

(5.7) T^"(A,k)= lim T(Aa^/^), 

n— >oo, a— 27r_R=na 

g«^(A,K)= lim a^"Q(Aa^/^). 

The eigenvalues of T^^{X, n), Q^^{\, n) are given by convergent infinite products due 
to (15.31) . In particular, it is easy to see from (15. 3p that the following asymptotics 
hold: 

(5.8) logg-(A,K) 27ri?--^ (-A^)^. 



sm „ 

21/ 
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Certainly, the limits exist for eigenvalues corresponding to any eigenvectors satisfy- 
ing ([53D. 

Now we turn to the operators Tm(C)'« + a — s), Qm(C)'^ + a — s) for which 
the eigenvectors of spin s are considered. The definitions of T'^'^{X,k, + a — s,s), 
Q'^'^{X, K + a — s,s) are the same as before. The important statement is that 

r^(A, K + a-s,s) = r^(A, k') , Q'%X, K + a-s,s) = Q''{X, k') , 

where k' is given by (14. 6p . The right hand sides are understood as "correct" analyt- 
ical continuations from spin sector. We have to explain two points: first, what is 
the reason that in the scaling limit the eigenvalues in the spin s sector equals analyt- 
ical continuations of those in the spin sector; second, what we mean by "correct" 
analytic continuations. The first point is simple: recall the discussion concerning 
the similarity of the equations (14. 5p . (14. 7p . The only difference between them was 
the number of Bethe roots, but this number is infinite in the scaling limit, so, the 
difference disappears. On the other hand the eigenvalue T'^'^(A, k) is a multi-valued 
function of k, so we have to explain the choice of its branch. At this point we refer to 
the semi-classical domain : u close to 1. We take a good branch at this domain, and 
then continue analytically. Notice that r(0, k) = 2cos(7rz/K). We require that intro- 
ducing s does not deviate us far from this value. This was the reason for choosing 
the definition (13. ip because with this definition we have 

r(0, K + a - s,s) = 2 cos(7rz/(K + a + 2^s)) , 

which stays close to 2cos(7rz/K) for all s if z/ is close to 1. 

From now on we shall often consider k and k' as arbitrary numbers implying the 
possibility of analytical continuation from values (14. 6p . 

Using Q^^{X, K,), Q^^{X, k') we obtain finite limits 

(5.9) p''''{X\k,, K,') = lim p{Xa''\hi,a, s), 

n— >oo, a— )-0, 27r/?=na 

(5.10) 4 oj^^{X, n' ,a) = lim ui^Xa" , fia'^\K,a, s) , 

n— >oo, 27r_R=na 

where k' is given by (14. 6p . and then for p^'^{X\k, k'), u^'^{X, p\k, K',a) the analytical 
continuation with respect to k' is used. 

According to Remark 14.11 (15.40 remains valid for any Bethe states in Matsubara 
direction. It has been already said that to Bethe states close to the ground states 
the scaling procedure applies. We shall argue later on that these vectors span the 
Verma module of chiral CFT. So, we would like to use the right hand side of (15. 4p 
in order to consider the scaling limit in the space direction 

ja = X finite. 

In our setting it amounts to considering the operators: 

(5.11) 2t*(A) = limt*(Aan, 

a-i>0 

2/3*(A) = limb*(Aa'^), 27*(A) = limc*(Aa'^) , 
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and 



$,(0) = limg2°^(°\ 

a->-0 



SO that (15 ■4p gives in the scahng hmit 

(5.12) ^r'{T*(A?) ■ ■ ■ T*{Xl)f3*{Xt) ■ ■ ■ /3*{Xt)Y{K) ■ ■ ■7*(Ar)($.(0))} 



fJp^^(A° I X det {co''i\t,\j\K,K',a)) 



i=l 



We understand the formula (15. 4 p as giving the expectation values of certain non- 
local operators making contact with quasi-local ones near ("^ = 1. After introducing 
a this point moves to = a~^'^, and in the scaling limit it goes further to = 
-|-oo. It should not be a surprise that this limit is described by CFT. To make this 
statement precise, we have to establish certain asymptotic properties of p^'^{X\k, k'), 
uj'^'^{X, fi\K, n',a) for A^,/i^ — > +oo. But first we shall need some information about 
the integrable structure of CFT. 

6. CFT ON A CYLINDER AND THREE POINT FUNCTIONS 

In this section we introduce our notation concerning CFT, and collect a few facts 
which will be relevant to the subsequent sections. 

Consider chiral CFT on a cylinder Cyl = C/27rzi?Z with circumference 2nR, the 
points X and x + 27riR being identified. Along with the local coordinate x, we shall 
also use the global coordinate 

X 

z = e R . 

The two points x = — oo, oo on the boundary of Cyl correspond respectively to the 
points 2; = oo, on the Riemann sphere. 
Let 



be the energy-momentum tensor in the coordinate x, where the /„'s satisfy the 
commutation relations of the Virasoro algebra with the central charge 

6z/2 

(6.1) c=' 



1-v 

The energy-momentum tensor in the coordinate 2, 



n-2 



is related to T(x) via the transformation rule T(2;)((i2;)^ = [T{x) — {c/12){z;x}){dx)'^. 
Here {z; x} denotes the Schwarzian derivative. In turn, T{x) is written as 



^ / \ 1 / X nx C 

Tix) = — > L„e- 

n=— 00 
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The Virasoro algebra acts on a local field 0{y) by the contour integral 

(6.2) {KO){y) = J ^Xx - yr^'T{x)0{y) , 

Cy 

where Cy encircles the point y anticlockwise. 

From now on, we fix a primary field (pAiy) with the scaling dimension A: 

(1o0a)(2/) = A(/.a(2/), MA)(y) = (n>0), 

and study the expectation values 

(6.3) (T(Xfc) ■ ■ ■r(xi)0A(l/))A+,A_ . 

The suffix indicates that we consider fl6.3p in the presence of two other primary fields 
inserted at x = ±oo. More precisely, we impose the boundary conditions 

is»^(-) = J? (^- - ^) 

inside the expectation values (16. Sp . where A± are the conformal dimensions of the 
inserted primary fields. For readers who prefer the language of representation theory, 
we are considering a highest weight vector |A+) at ^ = satisfying L„|A+) = 
5„,oA+|A+) {n > 0), and a co- vector (A_| at 2; = oo satisfying (A„|L„ = (5„,oA-(A_| 
(n<0). 

The singular part of fl6.3p is known from OPEs. In order to write them in the 
coordinate x, it is useful to introduce the function 

x(x) = ^coth(^)=f^^(|)" \ 

n=0 ^ ^ 

Here Bq = 1,-62 = 1/6,-64 = — 1/30,--- are the Bernoulli numbers. The main 
OPEs then read, as x ^ y, 

(6.5) r(x)T(y) = -^X"'{x -y)- ^^x'{x - y) + ^x{x - y) + 0(1) , 

(6.6) T(x)0a(z/) = -^^x'(x -y) + ^X{x -y) + 0(1) , 

where the prime stands for the derivative. 

The OPEs (16. 5p . (16. 6p . combined with (16. 4p and the behaviour 

X(x) = ±i + 0(e^^/^) (x^±oo). 
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allow us to write the conformal Ward-Takahashi identity which determines (16.3 
recursively: 

(T(Xfc) ■ ■ ■T{Xl)(f)Aiy))A+A- 
k 

^x"'ixi - Xj){T{xk) T(x2)0a(2/))a+,a_ 



c 



12R 



J=2 



+ (A+ - A^)^x{xi -y) + ^(A+ + A_) - ^}(T(xfc) ■ ■ •T(x2)0a(i/))a+,a_ 
From these one can extract, for example, 

((U0a)(i/))a„a_ _ J li^^^-"^- -12^- n{n-%Rn ^ 

{n : odd). 



(0a(i/))a+,a_ I (A+-A-)gn-i 

(n- l)!i?" 



In general, the normalised three point function of any particular descendant 

((!_„,, ■ ■■1-„i0a)(2/))a+,a_ 



(6.7) 



(0a(i/)) 



A+,A_ 



is determined as a polynomial in A, A+, A_. 

In writing these formulas, we have tacitly assumed that (0a(2/))a+,a_ is non- 
trivial. Actually, in the theory with c < 1, for a given generic value of A there 
is a discrete but an infinite collection of such A+, A_. In view of the polynomial 
dependence mentioned above, one can regard (16. 7p as a linear functional, defined for 
arbitrary A, A+, A_, on the Verma module consisting of the descendants of 0a (l/)- 
Henceforth we shall adopt this point of view. In later sections we shall use the 
parametrisation A = A^,, A+ = A^+i, A_ = A_k'+1; and write this functional as 



(6.8) Zr'{x{y)] 



{x{y)) 



Ak+i,A_ 



W(^/))a«+i,A_„,+i 



where X{y) is a descendant of (pAciv)- 

In [11], A. Zamolodchikov introduced the local integrals of motion which survive 
the 01,3-perturbation of CFT. They act on local operators as 

(6.9) (i,^_,0){y)= [ ^Mn{x)0{y) {n>l). 

Jcy 2vr2 

The densities h2n{y) are certain descendants of the identity operator /. The simplest 
examples are 

(6.10) h2{x) = {\-2l){x) = T(x), h,{x) = (I-2T)(X) , 
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for which we have 

oo 

(iiO)(y) = (l-iO)(i/), {hO){y) = 2j2ii^n-2\n-iO){y) . 

n=0 

In general, for a descendant h{x) of the identity operator, the three point function 

Cy 

is reducible to that of 0{y). Indeed, it can be rewritten as 

C{u) C{v) 

where C{u) {u G M) is a circle with the real part m, starting from u — irRi and 
ending at m + irRi. Choosing m -C Re y <^ w and using the boundary conditions at 
X — )■ ±00, one can show that each of them reduce to a constant. 
From the above remark it follows that 

(6.11) (i2n-l(0(2/)))A+,A_ = (/2V1 -4n-l) " (0(y))A^,A_ , 

where I^n-i denote the vacuum eigenvalues of the local integrals of motion on the 
Verma module with conformal dimension A±. Their explicit formulas for small n 
can be found in [9] (see also section [TOl f ll0.18p -( 110.20p below). Notice that, in the 
special case A+ = A_, the three point function vanishes on the image of the local 
integrals. 

As mentioned in the introduction, we accept the conjectural statement that the 
Verma module is spanned by the elements 

i2A;i-l ■ ■ ■ i2fcp-ll-2mi ' ' ' l-2m,(0a(O)) . 

Formula (16. lip tells that for the computation of the linear functional (16. 7p it suffices 
to consider the descendants by the even Virasoro generators {l_2n}n>i- 

Before closing this section, let us comment on a point which could be a source of 
confusion. The local integrals of motion arise in two different ways. At the boundary 
of the cylinder, they appear as the operators hn-i constructed from the modes L„ 
of the energy-momentum tensor in the coordinate e~^/^. These operators preserve 
the subspace of the Verma module of a given degree and can be diagonalised. In the 
classical limit, the eigenvalues of /2n-i correspond to the values of the integrals of 
motion on quasi-periodic solutions to the KdV hierarchy. In contrast, the action of 
the integrals of motion on local fields i2n-i are constructed from the modes 1„ in the 
coordinate x. Unlike in the first case, they do not commute with Iq (for example, 
the first integral of motion is l_i). They act as a creation part of the Heisenberg 
algebra. In the classical limit, they correspond to the action of the Hamiltonian 
vector fields generated by the local integrals of motion. So it does not make sense 
to talk about their diagonalisation. 
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7. Brief review of BLZ 

In a series of papers [9l [TOl HE] Bazhanov, Lukyanov and Zamolodchikov (BLZ) 
studied the integrable structure of the chiral CFT on a circle. We shall recall these 
results briefly, since they are quite relevant to us. 

It is convenient to write the energy-momentum tensor in terms of the chiral boson 

The operators P, Q are canonically conjugate, and the a„'s satisfy the Heisenberg 
algebra 

5 Vj — 2 ' V^rmO^n] — ^ Om+nfl ■ 

The energy-momentum tensor is expressed as 

(1 - u)T{x) =: ^\xf : +u^"{x) - \^ . 
The chiral vertex operator 

2 2 

(7.1) 0„(x) = e~3&t : eT^"^(^) : 

is a primary field of scaling dimension 

^ ° 4(1 -z/) ■ 

We note that the parameter ff^ used in [10] is identified as 

/3^ = 1-^, 

hence their g = e^'^^'^ is our — g^^. We have also changed the sign of P, and the 
normal ordering convention in [10] to 

2 2 
g X 

which results in the appearance of a scalar factor e « in (17. ip . 

The main object studied by BLZ is the universal monodromy matrix in CFT. 
It is an element of C^g(b^) ® Ah, where ?7g(b+) is the Borel subalgebra of \Jq{5\2) 
generated by eo,ei,/ii, and Ah is the algebra generated by P,Q,an (the suffix H 
stands for Heisenberg). Set 

2P 

(7.3) Xu{x) = eo®V+{x) + ei®V4x), :Kn = -- , 

1 — u 

(7.4) V±{x) = e-^^-^^Ti : e^^^^") : , 

so that [^Kh, Kt(3;)] = ±2V±(x). The universal monodromy matrix is defined to be 



(7.5) 



Th(A) = Texp a y %^{-%y)dy q 



2-kR 
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Here 7 exp stands for the path ordered exponential. Formula (17.51) is understood as 
a power series in A. The integrals in each term converge in the domain 1/2 < u < 1. 
Otherwise divergences occur and a regularisation is needed. 

We have considered two maps: Ug{b+) — t- End(V^) with two-dimensional Va, and 
t/g(b+) — !■ Osca with Osca being the g-oscillator algebra (see [2] for the notation). 
The images of 'Jh(A) under these maps are denoted by Ta^H(A) and Ta,h(A). Then 
following [IQ\ we define 



(7.6) T^^^(A) = Tr„ (T„,H(A)e-2-('^^®^)^ 

qCft^^) = A^(1 - e2-^^)Tr^ (T^,H(A)e^-^(^-^^)) 



There is a slight difference with \10\ due to different notation for the g-oscillator 
algebra. 

These operators satisfy the Baxter equation 

(7.7) Ti^^mTW = QS^"(Ag-^) + (A^) . 

An important property of these transfer matrices is that they commute with the 
local integrals of motion. The first local integral of motion is nothing but Lq — 
which commutes with the transfer matrices as mentioned above. Hence each of 
their eigenstate on a Verma module belongs to the subspace of a definite degree. 
In particular, the highest weight vector of the Verma module (primary field) is an 
eigenvector. 

Actually, the local integrals of motion are all encoded in the transfer matrix 
T^^^(A). The latter is known [9] to be an entire function of A^. One of the main 
statements of [9] is that it has the following asymptotics for A^ — )■ oo, A^ ^ I^<0; 

CXI 

(7.8) log(T^^^(A)) ~ RCo X^ + Y. CrX^^hn-i , 

n=l 



where C„ are known constants which can be found in [TO], they can be also extracted 
from section [10] of the present paper. For the moment the only point relevant to us 
is the fact that 

Ci < 0. 

This means that for sufficiently large A^ the highest weight vector is the eigenvector 
of 7h^"'"(A) with the maximal absolute value. The eigenvalue of Lq on the highest 
weight vector equals 

One may wonder why the asymptotic expansion ( 17. 8p is given as a series in the 
fractional power A^. As explained in [9], the reason is that A^ is a dimensionful 
quantity having the dimension of the inverse length. Indeed, consider the L-operator 
(17. 3p . It must have the dimension of the inverse length in order that the exponential 
in ( 17. 5p be dimensionless. But the operators V±{x) carry the anomalous dimension 
1 — V. So, clearly the dimension of A equals v. 
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We shall be interested in the Bethe roots, which are the zeros = Xf^ of <5h^^(A). 
They behave as 

Xl = 0{n^''), n^oo. 

Of equal significance are the zeros = /i^ of T^^^(A). The eigenvalue corresponding 
to the primary field has the characteristic property [10] that all A^ > > /i| for all 

8. CONFORMAL LIMIT IN THE MATSUBARA DIRECTION 

Let us return to the XXZ model. Using the notation introduced in section [5l we 
write the Baxter equation 

(8.1) TM(Aa'', HcjQMiXa", k) = a{Xa'^)QMiXa'' , n) + d{\a'')Qm{Xa!' k) , 
where 

a(Aa^) = (1 - ga^'^A^)", d{Xa^) = (1 - g^^a^^A^)" . 

We are interested in the maximal eigenvector \n) of T-^{\oy , k). 

We want to consider the limit n — )■ oo, a — )■ 0, while keeping na = 2t[R ■ C and A 
fixed. In this limit, for 1/2 < i/ < 1, 

a(A) ^ 1, rf(A)^l, 

so, if we identify 

(8.2) UK = -2P. 

the Baxter equation fl8.ip turns into (17. 7p . 

Now we want to fix the constant C in order to make the equivalence between Q^'^ 
and Q^^"'" exact. We had the asymptotics (15.81) . On the other hand, it is known [10] 
that 

logQ^^^(A, K)^R.±T (i^) r (l - i-) T{u)^ (-A^)^ . 



So, comparing we see that the agreement is exact if 

r (— 

(8.3) C= ^ ^".{. Tiu)^ . 

2v/^r(^) 

So, we come to 

(8.4) Q^^(A,«:) = g^^^(A)|^_^ . 

Let us argue that the vector \k) goes to the primary field with the dimension 

= 4(i4i;) - 1) ■ 

On the lattice the maximal eigenvector |k) is defined by the requirement that the 
eigenvalue T{1,k) is of maximal absolute value. In the scaling limit this corresponds 
to the requirement that T^^'^IX, k) is maximal for A^ large and positive. But the 
asymptotic behaviour of the BLZ transfer matrix is given by (17. 8p . so in the domain 

(8.5) -— <argA< — 
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the maximal eigenvalue corresponds to the primary field. Comparing with (17. 2p and 
(17. 9p , we find that in the picture of section [H] we have 

A+ = A,+i . 

Hence the boundary conditions at +oo is described by 0^+1 (+oo). Similarly, consid- 
ering the left Matsubara transfer matrix we find that its ground state corresponds 
to the scaling dimension 

(8.6) A_ = A_.,+i . 

In other words, in the picture of section [6] we have the left boundary condition 
described by the primary field cxd). 

The argument above is not rigourous, because it involves two limits which are 
a priori non-commutative. Nevertheless we believe it makes sense because of inte- 
grability, which stipulates that \k.) is an eigenvector of T^i^o,'^, i^) for all values of 
A. As a supporting argument, we quote from [T7] a knowledge that the eigenvalue 
T(C, k) of the lattice transfer matrix for r = g^/^ has maximal absolute value in the 
range |C| = 1, — 7rz//2 < arg^ < 7rz//2. This agrees exactly with (18. 5p . 

We call the previous reasoning a macroscopic one. For completeness let us give a 
less formal, microscopic derivation providing at the same time a constant which is 
important for physics. Consider the Matsubara transfer matrix T]v[(Aa'^, k), which 
is given explicitly by 



where 



.7) TM(Aa'^, k) = Tij (^Lj^ni^o.") ■ ■ ■ Lj^i(\a'')q 



m 



-(g — q ^)\a''a^ g2'^m _ a'^'^X'^q 
Let us make the gauge transformation 

(8.8) L,,^iX-an = g-^-^S-^-?4„.(Aa'^)gHE-T^-i 



3 



where 
2j,m(Aa'') 



1 - a^-A^g-- -(g - g-^) Aa-^g-^+E^il' a„ 

(g - q-')Xa''q--2-^'^=^' 1 - a^'^A^g'-- 

Now we recall known formulae concerning the continuous limit of the XXZ chain 
[T8] . A very accurate account of this matter is given in Lukyanov's paper [19]. 
Notice, however, that the Hamiltonian in |T9] differs from ours by a similarity trans- 
formation with the operator U = Ylm odd'^m- Having this in mind we rewrite the 
main order formulae for n ^ oo, y = ma (formulae (2.19) in [19]) as follows. 

(8.9) al ^ . . dy i^i-ty) - ^{-ty)) , 

(-1)™/f72 : e±^^^'"'^^^^"'^^^ : • 

Here ^p{x), V^(^) ^-^^ ^"wo chiral bosons with the same normalisation as in section [71 
The fractional power of a in the second formula is needed in order to compensate 
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the anomalous dimension of : e^'-'^^"*^^"'^^"*^^-' :, and F is related to the one point 
function of the latter [19]. From these formulae we see that 

(8.10) iy± = gTE-T^'^?^± ^ VZ a'-'^V^i-ty) , 

where V±{x) are the chiral vertex operators (17. 4p . The power of a changed due to 
a normal reordering, while the constant Z is obviously related to the asymptotical 
behaviour at m — )■ oo of the following two-point function for XXZ model: 

We do not know a direct way to fix this constant, however, our construction allows 
an indirect one. Indeed, in order to have complete agreement with CFT on this 
microscopic level we need that 

(8.11) 2i,m(Ar) = 1 + aXXj^ni-iy) + 0(a^^) , 
which would imply 

(8.12) Lj^n{\a'')---Lj^i{Xa'') ^ 7exp(^\ Xj^n{-iy)dy^ , 
where we have set 

X„h(x) = tq--^ {a+V^{x) + aTV+{x)) . 

So the microscopic picture agrees with the macroscopic one if 

1 

Z 



4sin2(7rz/)C2^ ' 
Altogether we obtain from (18.71) 

p2ttR 

(8.13) TM(Aa^/^) ^ Tr^- [e"*'^''Texp(A / %j^ui-iy)dy)], 

Jo 

giving rise to the BLZ transfer matrix (I7.6p with the identification (18.21) . In partic- 
ular, the second chirality decouples. 

9. CONFORMAL LIMIT IN THE SPACE DIRECTION 
Let us return to the formula 
(9.1) {t*(A?) ■ ■ ■ T*{XlW{Xt) ■ ■ ■ f3*{Xt)Y{K) ■ ■ ■ 7*(Ar) } 

= n P^^(A°|«:, k') det (u;-(A+, Xj\^, k\ a))^_^.^^ ...^^ , 

i=l 

We have seen that the functions in the right hand side are defined through the 
eigenvalues of the BLZ transfer matrix on the primary fields Thus 
the right hand side of (19.11) is defined. We want to interpret the left hand side of 
this equation. Our arguments are far from being mathematically rigourous, so, we 
formulate our statement conjecture. 
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Conjecture. Asymptotics of (19. ip for A^, A° — )■ cxd describes the expectation values 
of descendants for chiral CFT with c = 1 — of the primary field 0q,(O) inserted 
on the cylinder with the asymptotic conditions described by o-nd 0^+1 ■ 

Recall that we start with k, k', a satisfying (I4.6p . and then continue analytically. 
The three point function of the operators 0_k'+i(— C)o), 0q,(O), (^k+i(oo) does not 
vanish because (14.61) can be rewritten as 

{-k' + 1) + a + (/€ + 1) = 2 - 2^s , 

which coincides in our normalisation with the Dotsenko-Fateev condition [7] for one 
type of screening operators condition. We do not know if the second set of screening 
operators can be defined starting from the lattice model. 

In the present section we shall first present qualitative arguments in favour of 
this conjecture, and then explain how it can be verified quantitatively. The actual 
verification will be done in sections [TTl [121 

Consider the primary field g2°'^(°) on the lattice. Making the scaling limit in 
horizontal direction on the cylinder in the same way as it was done in the vertical 
one we conclude that this operator turns into $q,(0) = 0q,(O) ® 0q,(O). Typical 
operators in the space "\Na~s,s are of the form (f^°'~^^^'^^^cr^^- ■ ■ a^^O (s > 0) or 
^2(a-s)S(o)^- . . . ^-^Q ^ where is spinless. Then the same bosonisation 
formulae as (18. 9p 

imply that 

Scaling limit {'Wa-s,s) C V„ , ji^^ ® V_„ . 

So, the operators t*(A), /3*(A), 7* (A) do not change the Verma module for the 
second chiarlity. This is one reason to assume that they do not act on it at all. Let 
us give the first evidence for this claim. 

Consider the operator t*(A). It originates from its counterpart on the lattice, 
t*(C). We know that close to = 1 the operator t*{Q describes the adjoint action 
of XXZ local integrals of motion. So, we expect the same in the CFT. Using the 
BLZ formulae (17. 8p we see that 

oo 

(9.2) logp^^(A|/t,/t') ~ 5^A-'^C„(/2„-l(/€) -/2n-l(/t')) , 

n=l 

which implies together with (16. lip that 

oo 

r*(A) ~ exp(^^ A^^CJan-i) • 

n=l 

So, as it has been expected, the local operators are extracted from the action of 
T*(A) in the asymptotics at A — )• oo as coefficients of fractional degrees A^ which 
has the dimension of inverse length. Certainly, this exercise is quite tautological, 
and we would not write this paper if this were the only thing we can do. But it 
demonstrates the chiral nature of our operators. 
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In the appendix we prove the general statement that uj'^'^{X, ^\k,, K,',a) has the 
following asymptotics for A^,/i^ — )■ +00 

00 

(9.3) Ci;'"^(A, fj,\K,, k', a) ~ \/ ff^{X\K, k')^J p^'^{iJi\K^ k') ^ A ~ p ^~Uij{K, k.', a) . 

The proof will be given for the primary fields as asymptotical condi- 
tions, but it can be generalised to arbitrary descendants. So, in the weak sense the 
following operators are defined: 



(9.4) log(r*(A))^X^r;^._iA-^, 



2i-i 

^ 00 

(9.5) f3*{X) ~ y /3* _iA-^ 



00 

2j-l 



which act between different Verma modules: 



* 



Consider the Verma module V^. We obtain different elements of this module 
by operators T2j_i acting on the primary field , and by the same number of 
operators and 72j_i acting further. Due to the completeness in the lattice 

case in this way we obtain linearly independent vectors from Vq. Counting the 
characters we see that the entire Verma module is created in this way. Indeed, from 
a combinatorial point of view, r* is one odd boson, and /3*, 7* are Gross- Neveu 
fermions which in uncharged sector produce one even boson. 

No we shall proceed to computation of the coefficients Uij^K, k', a) and comparing 
them with the three-point functions of CFT. Since the operators t*(A) is already 
settled we shall consider the which means to ignore the image of the 

actions by the local integrals of motion i2n-i in the Verma module V^. Acting on 
the primary field 0^(0), the even generators of the Virasoro algebra, l_2fc, create the 
quotient space of the Verma module by these descendants. 

10. Asymptotics of logo*''' 
In this section we study the asymptotic behaviour of the function 

.c., ^ Q"%M, ^) 

as A^ — i- 00. Following closely the analysis developed in [10], we give a recursive 
algorithm for determining the coefficients of the asymptotic expansion. 
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It is known (see ^0\, (3.17) and (3.23)) that for large k the smallest Bethe root 
behaves as ~ c(z/)k^'', where 

(10.1) = (2]?) ' 

5 = log u — {1 — u) log(l — u) . 
The main technical idea in [10] is to consider the limit 



)2 

A^, K — J- 00 , keeping t = c{v)~^^^ fixed. 

Henceforth we change the variable from A to t and write 

F{t,K) = loga'^^(A,/€) . 

This function is to be determined from the non-linear integral equation. In order 
to write the equation, it is convenient to redefine some functions in terms of the 
variables t, u. We use 

K{t) = — \(t_l+l_t±l+l 
^ ' 2m 2 \tq^ - 1 tq-^ - 1^ 

We also use i?(t, u) to represent the following resolvent kernel. 

/di) 
—R{t, v)K{y/u) = K{t/u), (t, M > 1). 

1 

An explicit formula for i?(t, u) will be given below. 

The non- linear integral equation for F(t, k) (t > 1) reads 

00 

/du 
K(t/u)F(u,K,)— = -27nuK 
u 



1 

e"=-oo e~"=-oo 

K(t/M)log(l + e^("''^))^- j Alt/M)log(l + e-^("''^)^— 
1 1 
where e is a small positive number. From Appendix A, we see that 

F(t, k) = -F+(argt, K)\t\^ + 0{\t\'^^), < argt < tt, 

F{t,K)= F^{aigt,K)\t\^ +0{\t\-^), -7r<argt<0, 

where F± is positive. We seek for the solution of (110. 2p in an asymptotic series in 

00 



u 



5 



F(t,«:)^^«:-2"+iF„(t) 



n=0 



Consider first the leading coefficient -Fo(^)- For t > 1, from (9.2) follows 
(10.3) ((/ - K)Fo)it) = -2'Kiu, 
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where K denotes the integral operator on the interval [1, oo) 

oo 

du 



Kfit) = / Kit/u)fiu) 



u 



Equation (110. 3p can be solved by the standard Wiener-Hopf technique. 
Quite generally, for a function f(t) let 

oo oo 

-oo 

denote the Mellin transform and its inverse transform. For the solution of (110. 3p we 
shall need 

^ sinh(2. - l)vrfc ^ 
sinh Tik 

along with the Riemann-Hilbert factorisation 

i-k{k) = s{k)-^s{-k)-\ 

^ Til + {l-u)tk)T{l/2 + tuk) ^,Sk 
T{l + ik)^y2n{l - u) 

where 5 is defined in (110. ip . 

The function S{k) is holomorphic on the lower half plane Imk < l/2z/, and for 
— i- oo behaves as S{k) = 1 + 0{k~^). If we demand that 

Foit) = const, +o(^r^^ (t ^ oo) , 

then (110. 3p admits a unique solution given by 

(10.4) Fo{t)= I dW^S{l)j0—, (t>l) 

where 

1 

/ 



2v/2(l-z/) 

The right hand side of (9.4) gives a continuous function Foit) on the half line (0, oo) 
such that Foit) = for < t < 1. However, the function F^it) for t > 1 can be 
analytically continued to the sector | argt| < 2(1 — z/)7r, rewriting the equation (9.3): 

(10.5) Fo(t) = -27rzi/ + (irFo)(t), 

(10.6) (irFo)(t)= / dlfS{l)K{l) 



1(1 + 



2u' 



It is also possible to check directly the consistency of the formulas ( 110. 4p and (110. 5p . 
The difference of two integrals (110. 4p and (110.61) has the only pole in the upper half 
plane Im / + ^ > at / = 0, where we pick up the residue — 27riz/. 
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Now we turn to the higher order terms. Similarly as above, the Wiener-Hopf 
method allows us to find R{t, u) for t > 1, 

oo oo 



— oo-oo 

The analytic continuation is given by 

oo oo 

R{t,u) = K{t/u) + [ [ ^'^f'u'"'S{l)S{m)k{l)k{m) ~' 



27r 27r I + m - iO 



From this follows 



oo 



R{t,u)= [ —f'S{l)k{l)R{l,u) 
J 27r 



where 

oo 

dm , -i 



(10.7) R{l,u)= ! —u'^S{m) 

J 27r 



/ + m — zO 



-oo 



= u-''S{l)-^ + f —- ^ u'^^Sim^kim). 

The last line shows that R{1, e^) is analytic near x = 0. Equation (110. 2p can be 
converted into 

(10.8) F{t,K) = KFo{t) 

-(1 i?(t,w)log(l + e^("''^))^- I i?(t,n)log(l + e-^('^''^))^). 
1 1 

Motivated by the formula (110.51) . let us set 

oo 

(10.9) F{t,K) = KFo{t) + j dlf^S{l)k{l){<il{l,K) - K^Hoil)) , 

— oo 

where 'if{l, k) has an asymptotic expansion, 

oo . ~ 

(10.10) *(/,«:) c^^«:-2"+iv^„(0, *o(/) 



+ t) 

We show below that each coefficient \E'„(/) {n > 1) can be determined as a polynomial 
in / by a purely algebraic procedure. 
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With a change of integration variable, (110. 8p is brought further into the form 



-ico+e 



(10.11) ^{l,^)-K^,{l) = -j-\^ I gi?(/,e^-/^'^)log(l + e^(^'^'^'=''^)) 



ioo+e 

+ [ ^i?(/,e— /^'^)log(l + e^^(^'^^''^''')^ 
J 27r 



Since log(l + e^^^"'*^-*) decays exponentially for ±lmu > 0, the asymptotics of 
the right hand side of (110.111) is completely determined from the behaviour of the 
integrand at x = 0. 

In order to develop a systematic expansion, let us first make a general remark. 
Consider a Fourier integral 



/oo 
-oo 



We assume that g{k) is the boundary value of a holomorphic function on the lower 
half plane Im A; < 0, satisfying the asymptotic expansion 

oo 

5((fc) ~ ^ Qniiky^ {k^oo, lmk<0). 

n=— no 

Then integration by parts shows that for any > we have 

N 

G{x)= g^27r6^-\x) + 27rJ2—^xl-' + R^, 

n=—no n=l 

where Rn = 0{x^) as x — )■ 0. Suppose further that G{x) can be prolonged analyt- 
ically around a; = 0. In this situation its Taylor expansion can be computed from 
the right hand side, discarding the delta function terms. The result is summarised 
in a compact form 

fOO 



/oo 
e'^''g{k)dk = 27riresfc[e^'=^^(fc)] 
-oo 



where resfc[- ■ ■] signifies the coefficient of A; ^ in the expansion at A; = oo. 

The above consideration applies to (I10.7p . and we obtain the Taylor expansion at 

X = 0, 

(10.12) R{1, e'^/^^) = res, [^j^S{h) . 

For the factor log(l + e^^"''^^), we proceed as follows. Set 

F(e^^/^^ k) = -271 [x - F{x, K,)) . 
Similarly as above, the Taylor expansion of F{x, /t) at x = is calculated as 

(10.13) F(x, k) = X + res,, [e-'^^/^"5(/i)2^(/i, n) 
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Actually the term x is cancelled by a term coming from K\l/o(/i), so that F{x, k) 
0(k~^). We rewrite the corresponding piece of the integrand as 



(10,14) ,og(l + e-<"-"-.-)) =5:^(^(4)"log(l + 

n=0 

Substituting (I1U.14I) . ffTirT^ into f lTinT]) . we arrive at 
(10.15) 

k) - iK^Q{l) 

oo 

V— ^ 1 f (ir ( yp-hx/fK. -, _ / (9\"1 



-2'JTX\ 



Here {■ ■ ■ }even (resp. {■ ■ • }odd) means the even (resp. odd) part in x. To evaluate 
the integral in (110.151) we need only to develop the integrand into a Taylor series 
and apply the formula 

oo 


In summary, the asymptotic expansion (110. lOp can be calculated order by order 
in from the set of equations (110. 15p and (110. 13p . 
The first few terms of the expansion read 

^ V 2^y V ~ ' 7z/(i-z/) ) IfKf ^ ■ ■ ■ 

In general, the coefficients have the structure 

(10.16) ^„(/) = n 6 - X ( Polynomial in / of degree n - 1) . 

From the knowledge of log a'^'^(A, k), it is straightforward to extract the asymptotic 
expansion of logT'^'^(A, k) [lOj : 

V2^ JR~^^^o T{1 - t{l - iy)l) V 



ill. — "! 

Cn/2n-l(/t)A ^ 



2i/ 

OO 

/ \ N 2n-l 

n=0 
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where 

;io.i7) = -^3t;7^^i^i^7t:^^(i - n^r 



~^r(l + ^(2n- 1))' 



and we set by definition I_i = R. Notice that Cq > while C„ < for n >1. 

The factors in fll0.16p ensure that at these special values of / the asymptotic series 
fllO.lOl) truncates. This has to be the case, because according to [lOj hn-iii^) > 1) 
are the vacuum eigenvalues of the integrals of motion which are polynomials in c 
and A^+i. 

For instance 

(10.18) /,(«:) = i(^A,+i-^), 

(10.19) ^3(.)-^/.(.)^-^/.(K) + ^, 

(10.20) = ih{^)h{^) - -^M^) + ^M^) - ^ 



^ ov y XV ; ■^v ^ 3507^4 ^ 181440/?5 ' 

We have verified upto n = 4 that (110.171) matches perfectly the formulas for l2n-i 
given in [9]. 

11. ASYMPTOTICS OF CO FOR K = k' 

In this section, we restrict our consideration to the case k = k', so that p^'^{X\k, k') = 
1. Our goal is to give an algorithm for deriving the asymptotic expansion of the 
function u^'^^X, k, a). 

We start from the representation 

(11.1) (^""""{X, K, a) = ^/left * /right + /left * -Rdress * /right) (A, /i) + Uo{X, fi\a) , 

where 

/ieft(A,/i,a) = ^6^4j{X/n,a), /right(A, yu, a) = 5~V^(A//i, a) , 

wo(A,/i|a) = 5^5~A~V(A//i,a) , 
and -Rdress denotes the resolvent for the integral equation 

(11.2) /?dress -Rdress ^ -^a • 

Here we have set 




/(A)^7(A)rfm(A) 



2 



dX 



2 



A2(l + a-(A,ft:)) ' 

and cr^ is a point lying between the smallest Bethe root and the largest zero of 
T'"^(A, k). Strictly speaking, ic and dm{X) have slightly different meaning than those 
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used in section [3l Since we use them here only locally, there should not be a fear of 
confusion. 

From now untill (lll.sp below, we shall work with the variables 



t = c{yY^\'lK'\ u = c(z/)-V7/t'' , 

and write 



27rz 2 V' " ' tq^ -I ^ ' ' tq~^ - 1 

Again we start with the leading order approximation as k — t- oo, where flll.2p 
becomes 

R(t,u,a)— / — R(v,u,a)K(t/v,a) = K(t/u,a) . 
Ji V 

This can be solved in the same manner as before, using 

^ ' sinh7r(A; + f) 

The only point worth noting is that in writing the Riemann-Hilbert factorisation 

1 - k{k, a) = S{k, a)-^S{-k, 2 - a)-\ 

T(l + (l-u)ik-^)T(^ + iuk) .„ 
^ ^ r(l + zfc-f)v^(l-z/)(i~°)/2 

we are naturally led to assume that 

< a < 2. 

So the naive symmetry (— fc, —a) — {k, a) of K{k, a) is replaced by the symmetry 
{k, a) — )■ {—k, 2 — a). With our normalisation of a, the reflection a —i- 2 — a is the 
usual one for CFT with c < 1. For the resolvent kernel we obtain the representation 

R(t, u, a) = K{t/u, a) 

poo pea 7, , ^ ^ _ • 

+ / / ——f'u'"'S{l,a)S{m,2 - a)K{l,a)K{m,2 - a)- . 

J-oJ-oc 2n 2tt l + m-iO 

The 'dressed' resolvent kernel Rdress(t,u) satisfies 

Rdress{t,u) - R{t,U,a) 



/ :; ?7 — TRit,v,a)Rdmssiv,u) / — ^Rit, 



V, a)Rdvessiv 
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Setting 
(11.3) 



RdTess{t,u) = K{t/u,a) 



+ 



°° dl dm 



f a)S{m, 2 - a)K{l, a)K{m, 2 - a)e{l 



27r 27r 



;ii.4) 

oo 

9(/, m|/t, a) ~ 0n(^5 m|a;)K~^" , 0o(^) ^) 



n=0 



—I 



I + m 



and repeating the analysis of the previous section, we arrive at the hnear recursion 
relation for the 0„(Z,m|a): 



2 ^ 1 



6(/, ml/t, a) — 0o(/, m) ^ — — / dx 

■' n=0 



res;/ 



-I'x/fK 



-S{l',2-a) 



X reSr 



-■m'x/fKQf^l 



S{m' ,a)Q{m' ,m\K,a) F{x,k,) 



X 



1 



dx/ J odd 1 + e 



The coefficients of the series f lll.4p can be calculated by Taylor expanding the inte- 
grand and applying 



d 



dx) l + e^'^^' 

The first non-trivial term reads 

1 i 1 

^ ^ ' l + m 2Av {fny 



mid - 2 



-m+n' 



({m — n + 1) 
('27r)'"-"+i 



Returning to the original variables A and /x, formula for u^'^{X, fi\K, K,a) can be 
obtained from flll.ip . We have 



2m 



11.5) 



X 



dldmS{l, a)S{m, 2 — q;)0(/ + iO, m|/€, a) 



S{k, a) 



T{-ik + ^)T{^ + iiyk) 



r(-?(l - u)k + 1)72^(1 - z/)(i-")/2 ' 



where l + iO is important only in Oo(^, fri). Notice that originally in R{t, u, a) we had 
rather / — iO. The change appeared due to addition of uo{\, fi) which explains the 
importance of this term. Picking the residues at the poles in the upper half plane, 



36 H. BOOS, M. JIMBO, T. MIWA AND F. SMIRNOV 

its asymptotics as t, -u — > oo can be calculated: 

(11.6) a;^^(A,/i|/€,/€,a) - V — -D2r_i(«)/^2s-i(2 - a) 

^-^ r + s — 1 

r,s=l 



2r-l 23-1 

X A - /i - i22r-i,2s-i(ft:, a) 



where 



,n.T) Z,.„_,(a.) ^ ' m-il - . , \„ 

(ri-l)!rf|+ (k^(2n-l)j 



and 

2r+2s-2 



i(2r-l) i(2s-l) 



2z/ 2z/ 



r + s — 1 / V2 fKV 

K, a I X 



z/ \ R 



The counterpart of the factorisation 01O.16P for \l'n(0 is the vanishing property 

'i(2r-l) i(2s-l) 



2v 2v 



a = (n > r + s). 



This ensures that the coefficients ^^2r-i,2s-i('^, are polynomials in A^+i, a and c. 
For instance, 



rii,3 (k, a) = — Jsf/t) H ^ H rAa =F -do 

3,1^ ' ^ '^^ ^ 6/?3 1^ ^ 144i?4 ^ 1080i?4 "^ 360i?4 

1,5 («:,«) = -hU) - + ( + ^^A„^ hU) 

5,1^ ' ^ R ^ AR^ ' \48R^ 360R^ J ^ ^ 

A3 13(c + 35),2 2c2 + 21c + 70, 
Al ^ — A„ 



1728i?6 90720i?6 " 60480/^6 
^ V 120^ ^^''^ ~ 1440^ ~ 7560^ " ' 



f^3,3(«:, «) = -/5(/.) - ^/3(-) +[^. + + Y440^ J h{K) 

1 5c-14 10c2 + 37c + 70, l/2c2 + c 
K - r.. K ...... ^« — Aa 



1728i?6 " 18144i?6 " 362880i?6 ° 36288i?6 

Here hn-iii^) are given in fll0.18p -f lTo:20|) . and 

z/(z/-2) 
I/- 1 



(11.8) da= \^ / (a-l) = |V(25-c)(24A» + l-c). 



These structures exhibit a remarkable consistency with our fermionic picture. 
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12. Final results and conclusions 

Now we clearly see the structure of our fermions in the CFT limit. They naturally 
split into two parts: 

(12.1) = D^^_,{a)PZ% llra-l = D2m-l{2 - «)7^^!! , 

the multipliers D2m-i{oi), -D2m-i(2 — a) absorb all the transcendental dependence 
on a, the operators /32m-i! l2m-i purely CFT-objects. 

The fermions act between different Verma modules. In order to stay in one Verma 
module it is convenient to introduce the bilinear combinations of fermions 

/even _/ , -, \ ^ ( aC¥T* CFT* , ^CFT* CFT*\ 

H^2m-l,2n-l — [m -\- U l)^ {P2m-ll2n-l P2n-1 / 2m-l J 

^2m-l,2n-l = '^a ^ ('^ + ^ ~ 1) 2 {(^2n-l T2m-1 ~ f^2m-1^2n-l ) • 

The Verma module has a basis consisting of the vectors 

(12.2) 12/01-1 • ■ ■ hkp-l^-2h , • ■ ■ ^-2lq {<Pa) ■ 

Conjecturally the same space is also created by the action of the i2fc-i's and fermions: 

(12.3) i2A,i-l ■ ■ ' hkp~l4^2mi~l,2ni-l ' ' ' 4^2mr-l,2nr-1^2mi~l,2ni-l4^2rns-l,2ns-l{^a) ■ 

For small degrees, the transition coefficients between fll2.2p and (112. 3p . modulo 
descendants of the i2A:-i, can be determined by taking the expectation values with 
K = k' and equating like powers of k. Abbreviating (pa and writing as A, we 
find 
(12.4) 

jeven — l2 , ^'^ ~ 1 
^^1,3 = ^-2 H g 

J odd 2 I 

<Pl,3 = 3 1-4, 

_ 3 C + 2-20A + 2CA 

"^''^ = -2 + 3(KT2) 

-5600A + 428cA - Gc^A + 2352A2 - 300cA2 + 12c2A2 + 896A3 - 32cA3 
^ 60A(A + 2) 

,odd_ 2A 56 - 52A-2c + 4cA 

^^•^ = aT2 + 5(aT2) ' 

^^'^ 6(A + 2) 



+ 



-6544A + 498cA - Sc^ A + 2152A2 - 314cA2 + lOc^A^ - 448A3 + 16cA=^ 



60A(A + 2) 



1-6 , 



At the next degree, there are 5 Virasoro descendants Ii2, l-4l?-2, I-4, I-6I-2, 
1_8, which are polynomials in A^+i of degree 4,2,1,1,0, respectively. With the 
data at hand, obtained from the primary field (p^+i, there remains one parameter 
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undetermined. This can be fixed considering the first descendent L_i0k+i which we 
hope to do in future. Nevertheless we have checked that the determinant, 

/qCFT* aCFT* CFT*_ CFT* 
Pi Ps /3 ll 1 

after subtracting a suitable multiple of lig, has the correct degree 2 in A^+i. We 
regard it as a further supporting evidence in favour of the fermionic structure. 
Let us pass to conclusions. 

We believe that the fermionic description will provide new results for the theory of 
integrable models. For example, there is an obvious similarity between our fermions 
and those introduced in [20] . With the formulae fll2.4p at hand, it should be possible 
to upgrade the qualitative description of form factors of descendants in [20j to a 
quantitative level. We hope to explain this in future works. Here, however, we 
would like to emphasise that, even for CFT, the fermionic description must give 
something completely new. Let us explain that. 

Consider the functional Z'^'^ with k = k'. It describes the three point function for 
descendants of (pa and two primary fields (/>k+i of equal dimension A_k+i = 

Ak+1. It was said several times that the construction generalises if we replace the 
asymptotic states described by 0^+1 , by any other eigenstate of the integrals of 
motion l2n-i- The only change is that the function uj is to be computed for the new 
asymptotic condition. It is assumed [HIIID] that the joint spectrum of l2n-i is simple, 
so, in this way we compute all the three-point functions for a descendant of 0q, and 
descendants of provided the latter are eigenstates of the integrals of 

motion. Notice that the descendant of (pn+i can be very deep in the Verma module. 
In that case the usual CFT computation is rather hard to perform. Let us be more 
precise appealing to the classical limit. 

In the classical limit z/ — )■ 1, the eigenstates of /2n-i are in correspondence with 
the periodic solutions of the KdV equation. Let us give some explanation about this 
point. 

Consider the classical KdV hierarchy with the second Poisson structure: 
{u{yi),u{y2)} = 2{u{yi) + u{y2))6\yi - 1/2) + ^'"{yi - ^2) • 

The integrability of the KdV equation is due to existence of the auxiliary linear 
problem: 

(5,2+w(i/))V^(i/,A) = AV(i/,A). 

We consider the periodic case u{y + 2ttR) = u{y). In this case one defines the 
monodromy matrix M (A) for the auxiliary linear problem in a standard way. Then 
the local integrals of motion in involution are found in the asymptotical expansion 
of T'=^(A) = TrM(a) for A^ ^ +00: 

00 

log(T^'(A)) ^ 27iR\ + ^C^ij2=i_^A-(2"-i) , 

n=l 

where 

= -VTT j . 
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The integrals of motion I2I1-1 are well-known functional of u{y). Here we use for 
them the normalisation of [5]. So, the classical limit is 

T{-iy) Y^^^^^ 

It brings the Virasoro commutation relations to the second Poisson structure of the 
KdV hierarchy, and ensures the finite limits 

(1 - uri,^^^ ^ . 

It is well known that periodic solutions of KdV are in correspondence with hyper- 
elliptic Riemann surfaces which are two-fold covering of the Riemann sphere of A^. 
In particular, the solution corresponding after the quantisation to the primary field 
corresponds to the Riemann surface of genus 0: 

2 \ 2 2 

From the point of view of classical theory, this is a completely trivial case which 
describes a constant solution of KdV. This case becomes non-trivial after the quan- 
tisation, because KdV is a theory with infinitely many degrees of freedom, and 
quantising the simplest classical solution one has to take into account infinitely 
many zero oscillations (see [21] for a relevant discussion). Still it is rather unpleas- 
ant to be able to quantise only trivial classical solutions. The consideration of usual, 
low-lying descendants of does not change the situation seriously: they describe 
excitations for the same classical solution. What are really interesting solutions in 
the classical case? They correspond to other Riemann surfaces. The simplest one is 
described by the elliptic curve: 

/.2 = (A^-A?)(A2-A^)(A^-A^). 

At the quantum level, this solution corresponds to the following distribution of the 
Bethe roots over the real axis in the plane of A^. Going from A^ = —00 we first 
have no Bethe roots. Then there is a large interval where the Bethe roots are dense. 
Then there is a large interval without the Bethe roots, wherein we find one or several 
zeros of T^'^{X, k). Then starting from certain point and up to A^ = 00, the Bethe 
roots are again dense. 

For a reader who is not familiar with periodic solutions of KdV, it is useful to 
think about this solution as a periodic analogue of one-soliton solution which we 
really obtain in the limit i? — )• 00. Everybody would agree that quantising only the 
trivial solutions when there are solitons around is a waste of possibility. 

Our fermionic construction gives a possibility to treat this kind of asymptotic 
states. Moreover, we suppose that the function u has a clear algebra-geometric 
meaning in the classical limit. We hope to return to all that in one of our future 
publications. 

Appendix A. General results on the asymptotics of co'"'{\, K',a) 

In this section we derive the asymptotic behaviour ( 19. 3 p of (X'^"^(A, k,', a) when 
A^, /i^ — i- 00. The main point of the argument is that in a certain domain, which we 



40 H. BOOS, M. JIMBO, T. MIWA AND F. SMIRNOV 

call A-domain, the expansion (19 .2^ of log p^^{X\k,, k') holds for both A and Ag~^, and 
by the cancellation due to 

A^ = -(Ag-i)^, 

we have p^^{X\k,, K')p^^{Xq^^\K,, k,') ~ 1. We shall suppress the arguments k, k' and 
a in p^^{X\k,k') and u^'^^X, p\k, K',a). We set 

P (^) = ^sn/x ' a (A) - 



T'^{X,k) ' ' ' g«^(Ag"i,K) ■ 

For u^'^{X,p), after simple computations we get: 

(A.l) W'''(A, p) = (/left -k /right + /left * -Rdress * /right) (A, /i) + 6^ 6^ A^^tp {X/ p, «) . 

The symbol t^t stands for integration over the contour 7 going clockwise around the 
zeros of Q'^'^{X, k) with the measure 

d9^ 



dm{e) 



Here again -k, dm{X) and -Rdress are slightly different than those used in section |3] or 
section [TTl but this should not cause any confusion. 

The measure dm{X) has simple poles at the zeros of Q^'^{X, k) and T^"^(A, k,'). For 
simplicity of presentation, we assume k and k' are close enough so that any zero of 
T^^{X, K,') is smaller than any zero of Q^^{X, k). 

In this section if we say /(A) ~ g{X) on some half line of arg(A^), it means 
/(A) = ^(A) + 0(|A|-^) for all there. 

From [TD] one concludes that 

loga^^(A) = -F+(arg(A2))|A|^ + OdAp^), < arg(A2) < vr, 
loga^^(A)= F_(arg(A2))|A|H0(|Ar^), -vr < arg(A2) < , 

where F±(arg(A^)) are some functions taking positive values in corresponding do- 
mains. Hence a'^^(A) decays rapidly in the upper half plane and grows rapidly in the 
lower half plane. 

Following [To] we write in the integral over the upper bank in flA.ip using 



1 + 0*^^(77) 1 + a'^=(?7) ' a'^'^ir]) ' 

in order to separate the rapidly decreasing part. To formalise the story we introduce 
the notation 

f'*'9 = f°9-f*9 

where 

7(1 

e "00 00 

/(AMA)^^^^, f*g = jf{X)g{X)dm{X). 
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Here is an arbitrary point lying between the largest zero of T(a, k') and the 
smallest zero of Q^^[\^ n), and the modified measure is 

rfA^ / 1 1 



^"^(^) A2p-(A) Vl + a-(Ae*o) ' 1 + a-(Ae-^o) 
Introduce the resolvent R by the equation 

R- Ko,o R = K^, 

and two "dressed" kernels: 

(A. 2) Fieft = /left + /left ° R , -^right = /right + Ro /right • 

where the functions /ieft(A,?7), /right(^?, A) are defined in f l3.10p . They are singular at 
r/^ = A^. According to our general prescription we understand real A^ in them as 
A^e~*° and then continue analytically. The equation for the resolvent takes the form 

-Rdress + -R * i?dress = R * -Rdress + ^dress = R i 

and the definition of u can be rewritten as 

UJ^^\X, fl) = (-Fieft * Fright + Ficft * -Rdress * -bright) (A, /i) , 

a;(2)(A,p) = (/i^ft o Fright) (A, /i) + 6^6;A^'^{X/fi,a) . 

Now we are ready to study the asymptotical behaviour. We shall consider A^ and 
//^ in the A-domain defined as follows : 7r(2i/ — 1) < arg(A^), arg(/i^) < vr. We 
prove the correct asymptotic behaviour there, then assume that it is valid for all 
—71 < arg(A^), arg(/i^) < tt. The latter assumption is not even necessary for our 
goals, but we do not see why it should not be true having in mind that the only 
infinite series of poles of u{X,fi) are the zeros of T^'^{X,k,) which accumulate to 
A2 = -oo. 

The importance of A-domain is due to the fact that in it 
(A.3) p(A)p(Ag-i)^l. 
Introduce the operation 

(A.4) Stf{X) = f{X) + p{X)f{Xq~'). 

Using the definitions it is not hard to show that for A^, /i^ in A-domain 

5+Fieft(A, T]) ~ 0, 5+Fright(r?, yu) ^ . 

These equations imply 

oo 

2k~l 



(A.5) Fieft(A,r/) ~ Vp(A)5Z^"^^i«ft. '^(^)' 

k=l 
oo 

Fright(r?,/^) ^ v^p(p)^/i"^ Fright, kiv) ■ 

k=l 

It is easy to argue that Fieft, k{v)y Fright, k{v) grow for — > oo as powers of t]. This is 
enough to ensure that the "connected part" uj^^\X, fi) has the desired asymptotics: 
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we just substitute the asymptotics flA.5|) into the formula for uj^^\X, fx) and observe 
that all the integrals converge because of exponential in rj decay of drh{7]). With 
the "disconnected part" u^'^\X, n) the situation is far more dehcate, studying it we 
shall understand the importance of the term 6^6'^A^'^ip{X/fi,a) in uj{X,fi). 

Let us evaluate the last term in uj^'^\X, fj.) considering A^,/i^ > cr^. Using the 
definition (12.101) it is easy to see that 

e*°oo 

J 2z/(l + (A/?7)^) ^TT^V 4i/(l - (A//i)^) 

The function Fright (?7, A*) allows analytical continuation with respect to r], so, we shall 
use it for all G M+. Substite /right = Fright — o Fright and compute the integral 

u 

We get after some straightforward computations 

a;(2)(A,;,) =a;(3)(A,/i)+a;W(A,/i), 

where 

(A.6) u^'\X,ii) = - [ \ , ■ Fright(r/,/x) '^"'^ 



2u{l + {Xhf.) 2mrt^' 

o 

(A.7) UJ^'\\^^) = -VP 1 6-,5i \ ■ F,,,Uv,f^h '^"'^ 

J 2z/ 1 - A/?7 - 





00 



'(l-(A^)^) "^"^^"'^^27rzr^V(^)' 

u ~ 

where the principal value refers to the pole at r/^ = /i^. 

For a;*^^^(A,/i) the asymptotics of the kind fl9.3p follows immediately from (lA.SP 
and 

'^A'^A^ ^ TT-O- 

'2z.(l + (A^)^) 
Consider ijj'^'^\X, jj). We check the equations 
(A.8) 5,V^)(A,/i)^0, 5>(^)(A,/i)^0. 

The first of them follows immediately from two facts. First, 

^tKK—. — ^ TV - . 

" 2z.(l - (A/r^)^) 

Second, writing explicitly 

(A.9) 616^^ ' Piv)-pW , l-pWpiv) 



2u{l - (A/ry)^) 2u{l - (A/r/)^) 2u{l + (A/r/)^) 
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and recalling the asymptotic expansion for p(A), p{ri), we see that asymptotically 
for — >• +00, rf — >• +00 the singularities in f lA.Qp disappear. Altogether we have 
for the asymptotics in both arguments 

(A.IO) — -TT ^ E A-'^r . 

2z/(l - (A/r/)^) , 



To prove the second equation in flA.SP it is not sufficient to use f lA.SP because 



Fright('7, /i) has simple poles at rf = /i^ and 7]"^ = [x^q^ which contribute to the 
analytic continuation — )■ pq~^. However, it is easy to see that the corresponding 
contributions cancel. 

Using the first of equations (lA.SP we get 



00 

a;(^)(A,/i)^V^5^A-^u;^^)(/i), 

m=l 

where due to (lA.lOp the functions O0m\l^) are given by convergent integrals. These 
functions satisfy 5'^Um\fi) ~ 0, and do not grow for /i^ — )■ +00. Hence a;*^^^(A,/i) 
has the asymptotics of the kind (19. 3p . 
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